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Abstract: Climate change is threatening natural ecosystems in the Earth, and arid regions of southern 
Africa are particularly exposed to further drying. Welwitschia mirabilis Hook. (Welwitschiaceae) is an unusual 
gymnosperm tree that is recognized as an icon of the Namib Desert, southern Africa. Many aspects of its 
biology were investigated in the past, with a special emphasis for its physiology and adaptations, but 
nothing is known about its potential sensitivity to current climate changes. In this study, we adopted an 
approach based on distribution data for W. mirabilis and ecological niche models for clarifying the 
species-climate interactions and for predicting the potential impacts of climate change on W. mirabilis 
populations in three well-separated sub-ranges (northern, southern and central) in northwestern Namibia, 
southern Africa. We evidenced that the populations occurring in the northern sub-range have peculiar 
climatic exigencies compared with those in the central and southern sub-ranges and are particularly 
exposed to the impact of climate change, which will consist of a substantial increase in temperature across 
the region. These impacts could be represented by demographic changes that should be detected and 
monitored detailedly to plan efficient measures for managing populations of this important species on the 
long-term scale. 
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1 Introduction 


Global climate change is recognized as one of the major threats to biodiversity and natural 
ecosystems. Extinction events, range shifts, and density variations at population and species 
levels have been widely documented in cold, temperate, and tropical ecosystems (Walther et al., 
2002; Parmesan and Yohe, 2003, Cristofari et al., 2018). In addition, further pressures on natural 
populations, driving to changes in species diversity and ecosystem functioning and service 
provision are expected for the years to come (Hole et al., 2009; Moritz and Agudo, 2013; Ge et al., 
2015). Afro-tropical species and ecosystems are expected to be particularly sensitive to climate 
change and to face substantial alterations (Thuiller et al., 2006a; Peterson et al., 2014; 
Blach-Overgaard et al., 2015; Midgley and Bond, 2015). For this reason, demographic studies 
have been carried out on some species in arid regions of southern Africa for detecting climate 
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change impacts (Foden et al., 2007; Guo et al., 2011). 

Ecological niche models have been extensively used to project the likely impacts of future 
climate change for different time slices. In very general terms, the ecological niche models are 
based on the definition of the species-environment interactions and their spatial explicit 
projection in the future climate (Guisan and Zimmermann, 2000). In recent years, ecological 
niche models underwent deep improvements both in the theory and methodology (Guisan and 
Thuiller, 2005). A large spectrum of modeling approaches and techniques have been developed 
for dealing with different data types and aims (Elith and Leathwick, 2009; Guisan et al., 2013). 
These approaches have been used for predicting species distribution in a wide range of biota and 
for targeting initiatives, allowing to improve conservation strategies (Bakkenes et al., 2002; 
Iverson and Prasad, 2002; Thuiller et al., 2004, 2006b). 

Arid regions of southern Africa host a high percentage of endemic species and a variety of 
remarkable taxa (Maggs et al., 1998). Welwitschia mirabilis Hook. is an endemic, unique species 
of the taxonomic family Welwitschiaceae, and also an unusual member of gymnosperms in the 
Namib Desert, southern Africa. The distribution of W. mirabilis is restricted to the central and 
northern Namib Desert, extending from the Kuiseb River in Namibia to the Nicolau River, north 
of Namibe, in Angola (Kers, 1967; Giess, 1969). It occurs in isolated populations, with highly 
variable density (Henschel and Seely, 2000). Although a great number of studies have dealt with 
the morphology, development, ecophysiology, ecology, and distribution of W. mirabilis, many 
important aspects of the biology of this iconic species still need further understanding (Henschel 
and Seely, 2000; Valverde et al., 2016; Kriiger et al., 2017). To the best of our knowledge, no 
assessment of the potential impacts of climate change on this important species is available to 
forecast the variation in environmental suitability in the future. 

The aim of this study is to fill this gap, predicting how the climatic conditions in the Namib 
Desert will satisfy environmental exigencies of W. mirabilis populations in the future. We adopted 
an approach based on ecological niche models for clarifying the climatic components of the 
ecological niche of W. mirabilis populations and for evaluating whether and where such climatic 
exigencies will be satisfied across its current geographic range. More specifically, we combined 
distribution data with current climate surfaces for exploring the species-climate relationship and 
for modeling the climate suitability of W mirabilis populations in the Namib Desert under the 
current and future conditions. Furthermore, we projected a defined relationship into a wide 
spectrum of climate change scenarios for estimating the extent to which the climatic suitability will 
change Namib Desert in the future. Our results are important to understand more about the ecology 
of this unusual plant (i.e., W. mirabilis), for identifying which populations could suffer the worst 
impacts of climate change, and for improving the management strategies of W. mirabilis fields. 


2 Study area and data collection 


Distribution data for W. mirabilis were downloaded from the 'Digital Atlas of Namibia’ (Bubenzer 
et al., 2004), which provides maps compiled by the "National Herbarium of Namibia’ based on the 
databases of the 'National Botanical Research Institute of Namibia’ and published in the 'Atlas of 
Namibia’ (Mendelsohn et al., 2002) for the ‘Ministry of Environment and Tourism’ (Fig. 1). These 
maps are provided as polygon shape files for a selection of six plant species that are important at 
the country scale. Polygons describe the species distributions in terms of their extent of 
occurrence (IUCN, 2012), encompassing all known localities of presence. We defined our study 
area as a bounding box, which is 2/3 larger than the latitudinal and longitudinal extent of the data 
polygons. 

The available data could describe W. mirabilis as spread in three well-separated sub-ranges in 
the Namib Desert, southern Africa (Fig. 1). The northern sub-range is a circular area in central 
Kaokoland at the border between Namib Desert, Nama Karoo, and savanna biomes. The central 
sub-range is a large area in Damaraland and Swakopmund districts laying over Namib Desert and 
Nama Karoo. The southern sub-range is an elongated area in the central Namib Desert, 
Swakopmund district, southwestern Africa. The three areas are roughly 100 km away from each 
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other and are separated by unoccupied zones. 
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Fig. 1 Distribution of Welwitschia mirabilis in the Namib Desert, southern Africa 


3 Materials and methods 


Climate data for the study area were downloaded from the 'WorldClim databank’ (version 1.4; 
www.worldclim.org) (Hijmans et al., 2005a), which consists of 19 bioclimatic variables with 
global coverage. These variables were derived from the monthly temperature and precipitation 
values, using measures for the period 1950-2000 from a great number of weather stations 
(>15,000) around the globe. Data are available at different resolutions (i.e., 10.0, 5.0, 2.5, and 0.5 
min of arc) as grid (raster) layers. In particular, we used climate surfaces with pixels of 2.5 
arc-min, corresponding to a resolution of approximately 5 km (Hijmans et al., 2005b). 

In order to describe the climatic niche of W. mirabilis populations in a way that is easy-to-read, 
we represented its geographic distribution in a two-dimensional environmental space. For doing 
this, we extracted values of annual mean temperature and mean annual precipitation for all the 
pixels contained in our study area and projected the extracted values in a scatterplot. Values 
corresponding to pixels contained in the extent of occurrence of W. mirabilis were evidenced by 
coloring the relative points according to their sub-range. This approach produced a visual 
representation of the climatic niche of the species in its individual sub-ranges in relation with the 
climatic availability of the study area. In addition, we utilized the same approach for visualizing 
the climatic shift expected in the three different sub-ranges according to the available scenarios. 

We adopted an approach based on ecological niche models to estimate the climatic suitability 
for W. mirabilis populations across the study area (Guisan and Zimmermann, 2000). Since W. 
mirabilis has different climatic niches in its different sub-ranges (see Results for details), we 
modeled the species distribution in its different sub-ranges separately. We generated 1000 random 
points across the study area and considered them as pseudo-presence points when they fell into 
the considered sub-range and as pseudo-absences otherwise. Climate variables were assigned to 
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the presence/absence points and selected for reducing multicollinearity by discarding those with 
variance inflation factor higher than 5 (Belsley et al., 1981). After this selection, seven variables 
were retained for modeling (i.e., bio3 (isothermality), bio4 (temperature seasonality), bio9 (mean 
temperature of driest quarter), biol4 (precipitation of driest month), biol5 (precipitation 
seasonality), biol7 (precipitation of driest quarter), and bio18 (precipitation of warmest quarter)). 
The individual projections for the different sub-ranges were merged by summing up the suitability 
values in each pixel. 

Because the choice of algorithm can influence the performance of ecological niche models, we 
adopted an ensemble forecasting approach in this study (Aratijo and New, 2007). We used eight 
different techniques and then chose the most common and performing methods, for modeling 
climatic suitability in the R package (R Development Core Team, 2016) Biomod2 (Thuiller et al., 
2016). In particular, we fitted Generalized Linear Model (GLM) (Dobson, 1990), Generalized 
Boosting Model (GBM) (Friedman et al., 2000), Classification Tree Analysis (CTA) (Breiman et 
al., 1984), Artificial Neural Network (ANN) (Ripley, 1996), Surface Range Envelop (SRE) 
(Busby, 1991), Flexible Discriminant Analysis (FDA) (Hastie et al., 1994), Multiple Adaptive 
Regression Splines (MARS) (Friedman, 1991), and Random Forest Model (RFM) (Breiman, 2001) 
on presence/absence data for each sub-range. Each model was validated by the ROC (receiver 
operating characteristic) value (Fielding and Bell, 1997), and a final consensus model was 
generated by calculating the weighted average of the individual models using the relative ROC 
values as weights. 

To predict the variation of climatic suitability of W. mirabilis populations under future 
conditions, we projected the fitted models into the available Global Climate Models (GCMs). 
These GCMs were elaborated by 19 research centres for the four RCPs (Representative 
Concentration Pathway) for 2050 (average for the period 2041—2060) and downscaled at the same 
spatial resolution as the current climatic surfaces (i.e., 2.5 arc-min) using WorldClim as baseline 
climate (Hijmans et al., 2005a). Overall, we considered 59 different GCMs and used the same 
model parameterizations as in the current conditions for projecting the climatic suitability in the 
future and for creating the final consensuses model from the ensemble of the individual models. 
As a result, we obtained one single spatial prediction of climatic suitability for W. mirabilis 
populations under future conditions. 


4 Results 


4.1 Climatic exigences 


The evident separation of the three sub-ranges has implications for the climatic niches of W. 
mirabilis populations. Overall, this species occurs in the areas with very low annual precipitation 
respect to the entire study area (Fig. 2). In particular, populations in the southern sub-range (blue 
in Fig. 2) occupy the driest zones, with a mean annual precipitation <50 mm and an annual mean 
temperature of 15°C—21°C. Populations in the central sub-range (green in Fig. 2) have a relatively 
large spectrum of conditions for mean annual precipitation, from very dry (close to 0 mm) to 
semi-arid (250-300 mm), and a preference for hot conditions (annual mean temperature between 
19°C and 22°C) with few exceptions. Populations in the northern sub-range (red in Fig. 2) 
encounter conditions of mean annual precipitation in the spectrum of the central sub-range 
(roughly between 100 and 200 mm) but evidently colder (annual mean temperature between 17°C 
and 20°C). 


4.2 Climatic suitability 


Species distribution models obtained very high validation scores, i.e., mean ROC values of 0.84 
(0.20), 0.95 (+0.03), and 0.85 (+0.14) for the northern, central, and southern sub-ranges, 
respectively. The prediction map under current climatic conditions (Fig. 3a) evidenced that the 
northern and central sub-ranges are divided by a wide unsuitable zone. Differently, the central and 
southern sub-ranges are connected in terms of climatic suitability. The prediction map under 
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Fig. 2 Visualization of the climatic niche of W. mirabilis populations in the three main sub-ranges as well as the 
whole study area 
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Fig.3 Climatic suitability for W. mirabilis populations under current (a) and future (b) conditions 


future conditions (Fig. 3b) shows a similar pattern, with central and southern sub-ranges 
connected and northern sub-range isolated. In addition, lower levels of climatic suitability are 
forecasted for the northern sub-range in the future. 

The variation of climatic suitability between current and future conditions is expected to occur 
especially along the borders of species distribution range (Fig. 4). This is particularly evident for 
the central and southern sub-ranges, where the forecasted suitability variation has positive values 
(green in Fig. 4) along the borders facing the centre of the continent and negative values (red in 
Fig. 4) along the borders facing the ocean. The northern sub-range makes an exception to this 
pattern. Indeed, the entire northern sub-range is expected to suffer a very strong reduction of 
climatic suitability for W. mirabilis populations. 


4.3 Climatic variation 


The variation of climatic suitability evidenced above corresponds to a substantial increase of 
temperature expected across the species distribution. According to the 59 available GCMs, the 
temperature is expected to become higher in the future than in the current in the three sub-ranges 
(Fig. 5). In particular, the available scenarios of future climate predict that the annual mean 
temperature will become almost completely different from current conditions in the northern 
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Fig. 4 Variation of climatic suitability for W. mirabilis populations between current and future conditions. Note 
that green shades represent suitability gains and red shades represent suitability losses. 


sub-range (Fig. 5a). Differently, in the central and southern sub-ranges (Figs. 5b and c, 
respectively), the climatic niche of W. mirabilis populations will be still satisfied, although the 
temperature will change. These evidences are in accordance with results from ecological niche 
models. 


5 Discussion 


As evidenced by our results, the W. mirabilis populations in the northern sub-range have a 
climatic niche substantially different from those in the central and southern sub-ranges. This is 
evident both by the visualization in the two-dimensional environmental space and by the outputs 
of the ecological niche models. Indeed, the northern sub-range is colder than the central sub-range 
and wetter than the southern sub-range at the same time (Fig. 2). In addition, a wide unsuitable 
area is present between northern and central sub-ranges (Fig. 3). On the contrary, the central and 
southern sub-ranges appear relatively similar, being represented by a unique, continuous cloud of 
points in Figure 2 and being geographically connected in terms of climatic suitability (Fig. 3). 
This ecological difference between northern and central-southern populations could suggest some 
level of isolation between these two units, someway in agreement with the proposal of two 
distinct taxa (Leuenberger, 2001). Reasons behind Leuenberger's (Leuenberger, 2001) proposal 
(i.e., different cone morphology) were confuted by a recent study of Jacobson et al. (2014). But, 
unfortunately, populations from the Namibian northern sub-range were not considered in this 
study. Similarly, a preliminary study of Welwitschia genetic variability also did not include 
populations from the Namibian northern sub-range or from Angola (Jacobson and Lester, 2003), 
leaving the issue of the intraspecific taxonomic arrangement open. 

W. mirabilis populations in the northern sub-range are more exposed to impacts of climate 
change than those in the central and southern sub-ranges. Our results show that climatic 
suitability under future conditions will remain roughly the same for the central and southern 
sub-ranges. Differently, the northern sub-range will become rather unsuitable in terms of 
temperature (Fig. 5), losing most of the current climatic suitability in the future (Fig. 3). This 
evidences the importance of including W. mirabilis populations in the northern sub-range into 
future studies, monitoring plans, and assessing campaigns. Unfortunately, the northern sub-range 
lays in a remote region, where good roads are rare and traveling time is expanded. As a 
consequence, this area was seldom included in activities of systematic data collection and 
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Fig. 5 Variation of climate conditions expected in the three sub-ranges of W. mirabilis. In the figure, the big 
circles are the current values of temperature and precipitation in the northern, central, and southern sub-ranges, 


while the small circles are the future values in the same areas according to the 59 available Global Climate 
Models. 


available information mainly regards species occurrence or abundance. In the light of our 
outcomes, demographic studies aimed to understand whether climate change is influencing the 
structures and dynamics of W. mirabilis populations should be carried out, especially in the 
northern sub-range. 

The loss of climatic suitability we evidenced does not mean that of W. mirabilis populations 
will go extinct in the northern sub-range. Individuals of this species are among the oldest living 
plants, with ages estimated up to 3000 years (Bornman et al., 1972). In addition, it seems that this 
taxon faced substantial climate changes in the past. Indeed, some traits of W. mirabilis are not 
typically associated with desert plants (e.g., C3 metabolism, high water demands), justifying the 
hypothesis that earlier habitats of W. mirabilis were more mesic than the current habitat in the 
Namib Desert (Jacobson and Lester, 2003). In addition, paleontologic evidences suggest that W. 
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mirabilis populations existed in the areas corresponding to southern Africa and western South 
America more than 105 million years ago, before the separation of the continents and the 
upheaval of the Great Escarpment (Jacobson and Lester, 2003). Since the upheaval of the Great 
Escarpment was associated to high rainfall in the region, it is plausible that the current fragmented 
distribution of W. mirabilis was due to aridification during the Tertiary and Quatenary (Jacobson 
and Lester, 2003). In this light, living W. mirabilis trees are probably perfectly able to cope with 
the reduction of climatic suitability. Differently, recruitment could be affected by climate change. 
Indeed, a certain amount of rain, in specific rainfall pattern, is necessary for seed germination and 
seedling development (Bornman et al., 1972). For this reason, females annually produce many 
seeds but recruitment events are rare (many years apart) driving to coetaneous groups of 
individuals (Henschel and Seely, 2000). Therefore, climate change could prevent the germination 
of new plants, influencing the population structure. 

Obviously, some aspects of our study could be improved with additional data and information. 
First, distribution data we used are limited to the Namibian range and exclude the northernmost 
area of presence in Angola. Unfortunately, information about the species distribution in Angola is 
very general and, to the best of our knowledge, neither precise localities of species presence nor 
polygons indicating the extent of occurrence of W. mirabilis has been published. Moreover, the 
Angolan portion of W. mirabilis range is relatively small, being restricted to the coastal area south 
of Namibe (Kers, 1967). Second, distribution data we used are not precise localities of presence 
but general extents of occurrence. Unfortunately, available information about precise sites 
occupied by W. mirabilis is not abundant enough for niche modeling. In addition, density and 
precision of such information are heterogeneous among sub-ranges. This would have introduced 
into ecological niche models an unpredictable amount of uncertainty. Furthermore, techniques for 
cross-scale modeling are available (Bombi and D'Amen, 2012), which allow to produce detailed 
suitability models from general distribution data (Bombi et al., 2012). Third, ecological niche 
models assume the focal species is in a status of static equilibrium with environment. This 
assumption introduces a certain level of uncertainty into the outputs when dealing with 
non-equilibrium situations (Pearman et al., 2008). Modelers are developing new methods for 
incorporating adaptation in the process (Bush et al., 2016) to improve the predictive performances 
in dynamic conditions. Nevertheless, these methods require additional information about the 
species adaptive potential that is unavailable for W. mirabilis. Fourth, we used the ensemble of all 
the available GCMs for projecting the species niche into the future conditions. Recent studies 
evidenced that this approach can underestimate the impacts of climate change on wild species 
(Bombi et al., 2017), including unrealistic models into the possible scenarios. Unfortunately, 
unavailable information on species responses would be necessary for reducing the uncertainty 
related to future climate predictions. Overall, further detailed data about species distribution and 
biology would be needed to overcome these issues, but the state of the art of our knowledge 
already allowed us to evidence an effect. 


6 Conclusions 


Welwitschia mirabilis is recognized as a crucial element of the Namib Desert ecosystems. This 
species has a complex role in the ecosystems where it occurs, representing an element of 
attraction for many species of vertebrates and arthropods. In addition, this plant species is 
considered as a symbol of the Namib Desert and contributes to the wellbeing of local human 
populations drawing tourists in the area. Therefore, the development of a management plan for W. 
mirabilis populations, especially for those in the northern sub-range, would be fundamental and 
this plan should take into account the fact that climate change could have potential impacts on this 
species. In this regard, further efforts for collecting additional data on the distribution, 
demography, and adaptive potential of W. mirabilis populations are required. 
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